Explanation of observed features of self-organization in traffic flow 
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Based on simulations with the "intelligent driver model" , a microscopic traffic model, we explain the 
recently discovered transition from free over "synchronized" traffic to stop-and-go patterns [B. S. 
Kerner, Phys. Rev. Lett. 81, 3797 (1998)]. We obtain a nearly quantitative agreement with 
empirical findings such as the "pinch effect", the flow-density diagram, the emergence of stop-and- 
go waves from nonhomogeneous congested traffic, and the dimensions of their wavelength. 
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During the last years, theoretical and empirical inves- 
tigations have identified different possible mechanisms 
for a phase transition from free traffic to stop-and-go 
traffic on freeways. This includes deterministic 0-^] 
and stochastic mechanisms |^|| as well as effects of in- 
homogeneities p|-p^|. In contrast, Kerner has recently 
described the detailled features of another transition 
to stop-and-go patterns [ fl3[ developing from "synchro- 
nized" congested traffic ]14| , |l5[ on German highways, 
which arc compatible with empirical findings on Dutch 
highways |j~6[ . 

In the following, we propose a quantitative explanation 
of these observations based on microsimulations with the 
"intelligent driver model" (IDM). In particular, we will 
show the possible coexistence of different traffic states 
along the road behind an inhomogeneity of traffic flow. 
It is, in upstream direction, associated with the sequence 
"homogeneous cong ested traffic" fn||ll| (which, in a mul- 
tilane model, is related to the observed synchronization 
among lanes |) 17 1) — * "inhomogeneous congested traf- 
fic" |ll| (corresponding to the so-called "pinch region" 
p3[ ) — > "stop-and-go traffic" , while we have free traffic 
flow downstream of the inhomogeneity. 

It will turn out that, in contrast to previously reported 
traffic phenomena, this phenomenon relies on the exis- 
tence of a sufficiently large density region of convectively 
stable traffic, in which traffic flow is unstable, but any 
perturbations are convected away in upstream direction. 
Furthermore, one needs a traffic model in which the re- 
sulting traffic flow inside of fully developed traffic jams 
is much lower (nearly zero) than in "synchronized" traf- 
fic. In particular, without suitable modifications (see be- 
low), this is not satisfied by the traffic model discussed in 
Ref. We also point out that, although the IDM has a 
unique flow-density relation in equilibrium, it reproduces 
the observed two-dimensional scattering of flow-density 
data at medium vehicle densities [ fbi]]l~3[l , even without 
assuming a mixture of different vehicle types |lq] . 

The IDM is a continuous, deterministic model, in 
which the acceleration of a vehicle a of length l a at po- 
sition x a (t) depends on its own velocity v a (t) as well as 



the gap s a (t) = [x a -i(t) — x a (t) — l a ] and the velocity 
difference Av a (t) = [v a (t)— v a -i(t)] to the vehicle (a— 1) 
in front: 



(1) 



According to this formula, the acceleration on a free road 
(meaning s a — > oo) is given by a[l — (v a /vo) s ], where a 
is the maximum acceleration and vq the desired velocity. 
The exponent S is typically between 1 and 5. It allows to 
describe that the realistic acceleration behavior of drivers 
lies between a constant acceleration a up to their desired 
velocity vq (5 — > oo) and an exponential acceleration be- 
havior (5=1). 

The braking term —a(s*/s a ) 2 depends Coulomb-like 
on the gap s a , as it is the case for the braking term 
of the microscopic Wiedemann model |19|. Therefore, 
the acceleration term is negligible, if the gap s a drops 
considerably below the "effective desired distance" s* . 
With the relation 



s*(v a , Av a ) = so + max ( v a T + ° , 
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(2) 



it is constructed in a way that drivers keep a minimum 
"jam distance" sq to a standing vehicle, plus an addi- 
tional safety distance v a T, where T is the safe time head- 
way in congested but moving traffic. 

The nonequilibrium term proportional to Av a 
reflects an "intelligent" braking strategy, according to 
which drivers restrict their deceleration to b in "nor- 
mal" situations (e.g., when approaching standing or 
slower vehicles from sufficiently large distances), but they 
brake harder when the situation becomes more criti- 
cal, i.e., when the anticipated "kinematic deceleration" 
(Aw) 2 / (2s a ), which is necessary to avoid a collision with 
a uniformly moving leading vehicle (v a -i = 0), exceeds 
b. Notice that the acceleration a is typically lower than 
the desired deceleration b, and that both acceleration pa- 
rameters do not influence the equilibrium flow-density re- 
lation ("fundamental diagram"). Since it turns out that 
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neither multilane effects nor different types of vehicles 
are relevant in the context of this study, we have as- 
sumed identical "driver-vehicle units" characterized by 
the realistic parameters vq = 120 km/h, 5 = 4, a = 0.6 
m/s 2 , b = 0.9 m/s 2 , so = 2 m, and T = 1.5 s, apart from 
a localized change of vo or T (see below) . For the vehicle 
length we use I = 5 m, but this value does not affect the 
dynamics. 

We simulated an open freeway section of 20 kilometer 
length for time intervals up to 120 minutes, of which we 
display the most interesting parts only. In addition, we 
assumed an inhomogeneity of traffic flow that will be re- 
sponsible for the transition from free to congested traffic, 
as described in Refs. |lo| , |i"T|| . However, as pointed out 
by Kerner |l3| , the self-organized patterns observed by 
him are not restricted to the vicinity of on-ramps. We 
will confirm this by simulating different kinds of inho- 
mogeneities (see Figs. [j] through ||) and comparing them 
with the injection of vehicles at on-ramps (see Fig. 

In Figure |l|, we have assumed an inhomogeneity cor- 
responding to a freeway section where people drive more 
carefully. This was modelled by setting the desired time 
headway from T = 1.5 s to T = 1.75 s between x — km 
and x = 0.3 km. In the simulations of Figs. || and ||[ we 
have reduced the desired velocity from vq = 120 km/h to 
vq = 80 km/h in the same region. As initial conditions, 
we assumed homogeneous free traffic in equilibrium at 
a flow of Qq = 1670 vehicles/h (Figs, [j] through ||) or 
1570 vehicles/h (Fig. [|). The actual initial conditions, 
however, are only relevant for a short time interval. At 
the upstream boundary, we assume that vehicles enter 
the freeway uniformly at a rate Qi n (t) = Qq + AQ(t) and 
drive with a velocity corresponding to free traffic in equi- 
librium. While in the simulation of Fig. [I], the breakdown 
of traffic flow is caused by exceeding the static freeway 
capacity at the inhomogeneity jy]], in Figs. || and || it 
is triggered by a triangularly shaped perturbation AQ(t) 
of the inflow @0 between t = 10 min and t = 20 
min with a maximum of 200 vehicles per hour and lane 
at t = 15 min. Furthermore, we used the "absorbing" 
downstream boundary condition v a = 0. To minimize 
simulation time, we integrated Eq. (Q) with a simple Eu- 
ler scheme using a coarse time discretization of At = 0.4 
s, and translated the vehicles in each step according to 
x a (t + At) = x a (t) +v a At+ ^v a (At) 2 . However, smaller 
values of At yielded nearly indistinguishable results. 

Figure ^ gives a representative overview of the simu- 
lation result by means of a spatiotemporal density plot. 
The resulting sequence of transitions is essentially the 
same as observed [jDj): After 10 minutes of free traffic, 
traffic breaks down near the inhomogeneity, resulting in 
homogeneous congested traffic at this location, that per- 
sists over a long time. Upstream of the inhomogeneity, 
small oscillations develop that travel further upstream 
and grow to stop-and-go waves of relatively short wave- 
lengths (about 0.8 km). Finally, these waves either dis- 



solve or merge to a few "wide jams" (in which traffic 
comes to a standstill) with typical distances of 2 km up 
to 5 km between them. Once the jams have formed, they 
persist and propagate upstream at a constant propaga- 
tion velocity without further changes of their shape. No 
new clusters develop between the jams. 

To compare our simulation results directly with the 
empirical data published by Kerner |l3| , we investigated 
the temporal evolution of the average velocity at six sub- 
sequent locations Dl through D6 that had the same dis- 
tances with respect to the inhomogeneity as in Ref. |l3] ] 
(see Fig. ^). The detector positions Dl through D4 are 
upstream of the inhomogeneity, D5 is directly at the in- 
homogeneity, and D6 is downstream of it. In contrast 
to the simulation of Fig. El the capacity drop at the 
inhomogeneity is so weak that free traffic is metastablc 
in the overall system. At the inhomogeneity (D5), one 
observes homogeneous congested ("synchronized") traf- 
fic, at D4 one sees small oscillations that, around D3, 
develop to stop-and-go waves of larger amplitude, and fi- 
nally to jams (at Dl and D2). In the downstream direc- 
tion, the congested traffic dissolves to free traffic (D6). 
Apart from irregularities in the measured data due to 
fluctuations, the curves are in (semi-)quantitative agree- 
ment with Kerner's empirical findings p3| . 

We also plot one-minute data of the six "detectors" 
in a flow-density diagram (Fig. ||), together with the 
equilibrium flow-density relation Q c (p) ("fundamental 
diagram"). The lower boundary of the data points at 
medium densities corresponds to the flow-density relation 
belonging to the downstream front of a fully developed 
jam. In agreement with observations, this line is the same 
for all jams and corresponds to a unique propagation ve- 
locity of their downstream fronts. The twodimensional 
region of points above this line relate to congested traffic 
at D3, D4, and D5. 

To understand this scenario, we need some basic re- 
sults about the stability of homogeneous traffic with re- 
spect to localized perturbations. Typically, there are four 
"critical" densities p C i with r>ci < p C 2 < Pc3 < PcA and 
the following properties (2^j2^|: For low and very high 
densities {p < p c \ or p > p C i), traffic is stable with re- 
spect to arbitrary perturbations, while for p C 2 < p < p C 3, 
it is linearly unstable. In the two density ranges in be- 
tween, homogeneous traffic is unstable only with respect 
to perturbations exceeding a certain critical amplitude 
Ap CI (p) ( "metastability" ) . Furthermore, there exists a 
range p cv < p < p c y, with p cv > p C 2, where traffic is 
linearly unstable, but convectively stable, i.e., all pertur- 
bations grow, but they are eventually convected away in 
upstream direction p2| . Actually, this range can be very 
large. For the IDM with the parameters used here, we 
have pcv = 50 vehicles/km and p C 3 = 100 vehicles/km. 

Let us assume that the inflow Qm has a value larger 
than Qci, and that a phase transition from free to con- 
gested traffic has occurred at the inhomogeneity. This 
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breakdown to congested traffic, which we identify with 
"synchronized traffic" , here [ ]io| , [iT||l^|]l| | , was explained 
in [Tc| ]. Since it is connected with a drop of the effec- 
tive freeway capacity to the self-organized outflow Qout 
from "synchronized" traffic (ST), the region of congested 
traffic will grow in upstream direction until the inflow 
Qin{t) — Q c {Pm{t)) to the freeway falls below the syn- 
chronized flow QgT in the congested region behind the 
inhomogeneity. Without an on-ramp flow Qramp per 
freeway lane, we have Qst = Qout, otherwise it is 
Qst = (<3out - Qramp) @- In_ contrast to the empir- 



ical findings explained in Ref. |10|, the synchronized flow 
must be so high, here, that it is linearly unstable, which 
implies Qst > Qe{pcz)- Therefore, small perturbations 
will grow to larger oscillations, which propagate in up- 
stream direction faster than the congested region grows. 
When the oscillations reach the metastable region of free 
traffic upstream of the inhomogeneity, oscillations with 
an amplitude below the critical amplitude Ap cr (/9j n ) will 
eventually disappear, the remaining ones will continue to 
grow until they are fully developed traffic jams. We point 
out that a sequence of such jams is sustained, because the 
propagation velocity v g and the outflow Q ut ~ Qc(Pci) 
from jams are characteristic constants p3| , p4| . Finally, if 
the inhomogeneity is such that the synchronized flow is 
convectively stable, i.e., Qst < Qcv = Qe(pcv), the per- 
turbations cannot propagate downstream (in contrast to 
small perturbations in free traffic). Hence, the front of 
dissolving traffic is smooth, then, and a small region of 
homogeneous congested traffic forms near the inhomo- 
geneity [Figs. @, @(a), and§. 

In Fig. ||, the flow-density relation of the downstream 
front of a fully developed jam (where the velocity is zero) 
corresponds to a straight line, the slope of which is v R [p3[ . 
Notice that this line (which in Fig. || and in Ref. |13[ is 
labelled by "J" ) lies considerably below the equilibrium 
curve Q e (p). On the other hand, traffic flow in the region 
of oscillating congested traffic is nearly in equilibrium 
as long as the oscillations are small. As the oscillations 
grow, the data points gradually approach the line J. This 
explains the observed twodimensionality of the congested 
part of the flow-density diagram. We point out that mix- 
tures of different types of driver-vehicle units lead to to 
further effects contributing to an even wider scattering 
of flow-density points in the congested regime [Q . 

In summary, we have shown that the emergence of 
stop-and-go waves out of synchronized traffic, their co- 
existence, and the twodimensional scattering of data in 
the congested part of the flow-density diagram can be 
explained in the framework of "standard" traffic mod- 
els that have a unique equilibrium flow-density relation. 
The necessary conditions are a metastability of traffic 
flow, a flow inside of traffic jams that is much lower than 
in synchronized congested traffic, and a sufficiently large 
density regime of linearly unstable traffic flow that is con- 



vectively stable. If "synchronized" traffic is linearly un- 
stable and free traffic upstream is metastable, upstream- 
moving perturbations will grow and, when their ampli- 
tudes became large enough, eventually form stop-and-go 
waves. To maintain this mechanism, the region of syn- 
chronized traffic behind the inhomogeneity must persist 
(i.e. it must not dissolve to stop-and-go waves), which 
is only the case if it is convectively stable. The twodi- 
mensionality of the congested branch of the flow-density 
diagram originates from the fact that the nonhomoge- 
neous congested states are not in equilibrium. 

The phenomenon should be widespread since it is trig- 
gered at relatively small inhomogeneities whenever traffic 
flow Qi n (t) exceeds a certain threshold Q c \. Note that it 
can be also simulated with other traffic models like the 
macroscopic gas- kinetic- based model 0| , if "frustration 
effects" are additionally taken into account (see Fig. |J) ■ 
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FIG. 3. Flow-density data (symbols) calculated from 
the microscopic simulation at the four detector locations dis- 
played in Fig. |^, and in addition at x = —3.7 km (detector 
D2) and at x — 1.2 km (D6). Also shown is the equilibrium 
flow-density relation (solid line), and the flow-density relation 
at the downstream fronts of fully developed jams (dashed). 
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FIG. 1. Spatiotemporal density plot illustrating the 
breakdown to "synchronized" traffic (smooth region of high 
density) near an inhomogeneity, and showing stop-and-go 
waves emanating from this region. Traffic flows in positive 
x-direction. The inhomogeneity corresponds to an increased 
safe time headway T between x — km and x = 0.3 km, 
reflecting more careful driving (see main text). Downstream 
of the inhomogeneity, vehicles accelerate into free traffic. 
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FIG. 2. Temporal evolution of the average velocity deter- 
mined from the individual velocities of the vehicles that pass 
cross sections of the freeway during one-minute intervals at 
the four "detector" positions Dl, D3, D4 (upstream of the 
inhomogeneity), and D5 (at the inhomogeneity). The naming 
of the detectors and their distances with respect to the inho- 
mogeneity are the same as in Ref . jl3| . The inhomogeneity is 
realized by a drop of the desired velocity (see the main text), 
(a) Breakdown to homogeneous synchronized traffic around 
t — 45 min. (b) Oscillating synchronized traffic ("pinch re- 
gion"), (c) developing stop-and-go waves, and (d) resulting 
traffic jams. 



FIG. 4. Spatiotemporal evolution of the traffic den- 
sity according to the gas-kinetic-based traffic model [ pl[ . 
The assumed inhomogeneity of traffic flow comes from an 
on-ramp of length 200 m with an inflow of 220 vehicles 
per hour and freeway lane. The inflow to the main road 
is 1570 vehicles per hour and lane, and the breakdown of 
traffic flow is triggered by a perturbation AQ(t) of the in- 
flow with a flow peak of 125 vehicles per hour and lane (see 
main text). The assumed model parameters are Vo = 120 
km/h, T — 1.5 s, r = 30 s, p max = 120 vehicles/km, and 
7 = 1.2, while the parameters for the variance prefactor J^] 
A(p) = A + AA{tanh[(p - p c )/Ap] + 1} are A = 0.008, 
AA = 0.02, p c = 0.27p max , and Ap = 0.1p m ax- In or- 
der to have a large region of linearly unstable but convec- 
tively stable traffic, we introduced a "resignation effect", i.e. 
a density-dependent reduction of the desired velocity Vq to 
Vo'O) = Vo - AV/{1 + eM(p'c - P)/Ap']} with AV = 0.9Vb, 
p' c = 0.45p ma x, and Ap = 0.1p max . 
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